GENERATION OF TWO-DIMENSIONAL WATER WAVES BY MOVING 

BOTTOM DISTURBANCES 
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Abstract. In this study we investigate the potential and limitations of the wave genera- 
tion by disturbances moving at the bottom. More precisely, we assume that the wavemaker 
is composed of an underwater object of a given shape which can be displaced according 
to a given trajectory. The practical question we address in this study is how to com- 
pute the wavemaker shape and its trajectory in order to generate a wave with prescribed 
characteristics? For the sake of simplicity we model the hydrodynamics by a generalized 
forced BBM equation. This practical problem is reformulated as a constrained nonlinear 
optimization problem. Some constraints are imposed in order to make practically feasible 
the computed solution. Finally, we show some numerical results to support our theoretical 
and algorithmic developments. 
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1. Introduction 

The problem of wave generation is complex and has many practical applications. On the 
scale of a laboratory wave tank a traditional wavemaker is composed of numerous paddles 
attached to a vertical wall and which can move independently according to some prescribed 
program. These systems have been successfully used to conduct laboratory experiments at 
least since late 60's [25, 37]. 

In this study we investigate theoretically and numerically the potential for practical 
applications of a different kind of wave making devices. Namely, the mechanism considered 
in this study is composed mainly of an underwater object which can be displaced along 
a portion of the bottom with the prescribed trajectory. In mathematical terms, we study 
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the wave excitation problem by moving forcing. Similar processes are known in physics 
under the name of autoresonance phenomena, thoroughly studied by L. Friedland and 
his collaborators [26, 27, 28]. 

Recently, this type of wavemakers have found an interesting application to the man-made 
surfing facilities [34]. This type of devices was proven to be successful to generate high 
quality waves for surfing far from the Oceans. Our main goal consists in providing the first 
elements of theoretical explanation of this process. The second objective of this study is 
to provide an efficient computational procedure to determine the underwater object shape 
and trajectory to generate a prescribed wave profile in a given portion of the wave tank. 

The problem of wave generation by moving bottom has been particularly studied in 
the context of tsunami wave generation. These extreme waves are caused by sea bed dis- 
placements due to an underwater earthquake [32, 9, 41, 19, 24] or a submarine landslide 
[60, 42, 61, 2]. However, the bottom motion produced by an underwater earthquake is 
mainly vertical, even if some effort has been made to include also the horizontal displace- 
ments components [55, 57, 56, 40, 23]. On the other hand, the wave making mechanism 
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studied here involves only the horizontal motion. Consequently, the knowledge of the 
tsunami wave community cannot be directly transposed to this problem. 

The wave propagation takes place in a shallow channel, so the long wave assumption 
can be adopted [59, 15]. However, weak dispersive and weak nonlinear effects should be 
included since the resulting wave observed in experiments has some common characteristics 
with a solitary wave. Consequently, as the base model we choose the classical Boussinesq 
system derived by D.H. Peregrine (1967) [43] and generalized later by T. Wu (1987) 
[62], who included the time-dependent bathymetry effects. In order to simplify further the 
problem, we assume the wave propagation to be unidirectional and, hence, we derive a 
generalized forced BBM equation [3]. This equation is then discretized with a high order 
finite volume method [4, 21, 12, 22]. Finally, the trajectory and the shape of the underwater 
object are optimized in order to minimize a cost-function under some practical feasibility 
constraints. 

From mathematical theoretical point of view our formulation can be seen as the con- 
trollability problem for the forced BBM equation [45]. Let us describe the main available 
results on the controllability of dispersive wave equations such as KdV [8, 36], BBM [3] 
and some Boussinesq- type systems [6]. 

The controllability of the KdV equation: 

Ut + u xxx + u x + uu x = 0, x e [0, L] ,t > 

is well studied in the literature. The controllability and stabilization properties were ob- 
tained by L. Russell & B.-Y. Zhang (1993) [50] for periodic boundary conditions with 
an internal control. The boundary control was investigated by the same authors later [51]. 
The controllability of the KdV equation with Dirichlet boundary conditions was studied 
in the following papers [46, 63, 44, 13, 48, 10, 30, 11, 31]. This list of references is not 
exhaustive. 

Let us briefly describe now some results on the controllability of the BBM equation. 
Rosier & Zhang (2012) [49] proved the Unique Continuation Property (UCP) for the 
BBM equation posed on a one dimensional torus T := M/27rZ: 

Ut ~ u txx + u x + uu x = 0, (1.1) 

i.e. for any open nonempty set u c T the only solution of (1.1) with 

u(x,t) =0 for (x,t) eux (0,T) 

is the trivial solution u = 0. Moreover, they proved the UCP also for BBM-type equations 
of the form 

u t - u txx + [f(u)] x = 0, 

where / 6 C 1 (1R),/(m) > for all u e R, and the only solution u e (S,5) of f(u) = is 
u = 0, for some number 5 > 0. Furthermore, they consider the following control problem: 

Ut - u txx + u x + uu x = a(x + ct)h(x, t), 

where a e C°° is given and h(x, t) is the control. They prove local exact controllability in 
H S (T) for any s > and global exact controllability in H S (T) for any s > 1. A necessary 
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Figure 2. Sketch of the physical domain with an underwater object moving 
along the bottom. 

and sufficient algebraic condition for approximate controllability of the BBM equation 
with homogeneous Dirichlet boundary conditions was given in [1]. The controllability of 
linearized BBM and KdV equations was studied in [47, 38, 65]. The controllability of 
a family of Boussinesq equations have been studied theoretically as well [64]. In [58], 
Touboul recently obtained controllability results for the heat and wave equations with a 
moving boundary. 

The present study is organized as follows. In Section 2 we derive the governing equation. 
Then, this model is analysed mathematically in Section 3. The results of some numerical 
simulations are presented in Section 4. Finally, in Section 5 we outline the main conclusions 
of this study. 



Consider an ideal incompressible fluid of constant density in a two-dimensional domain. 
The horizontal independent variable is denoted by x and the upward vertical one by y. The 
origin of the Cartesian coordinate system is chosen such that the line y = corresponds to 
the still water level. The fluid is bounded below by an impermeable bottom at y = -h(x,t) 
and above by an impermeable free surface at y = t](x,t). We assume that the total depth 
H(x,t) = h(x,t) + rj(x,t) remains positive H(x,t) ^ Hq > at all times t. The sketch 
of the physical domain is shown in Figure 2. The depth-averaged horizontal velocity is 
denoted by u(x, t) and the gravity acceleration by g. The fluid layer has the uniform depth 
d everywhere, which is perturbed only by a localized object, which can move along the 
bottom: 



2. Mathematical model 



h(x,t) =d-C(x,t), COM) =( (x-x Q (t)), 



(2.1) 



where the function Ca( x ) has a compact support and x = x (t) is the trajectory of its 
barycenter. The meaning of the segment X = [a, b] is explained in Section 3.1. 
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In 1987 T. Wu [62] derived the following Boussinesq-type system to study the generation 
of solitary waves by moving disturbances: 

Vt + ((h + V )u) x = -h t , (2.2) 

u t + uu x +gr] x = ^h(h t + (hu) x ) xt -^h 2 u xxt . (2.3) 

This system represents a further generalization of the classical Boussinesq equations derived 
by D.H. Peregrine (1967) [43] for the case of the moving bottom h(x,t). In our work we 
take the system (2.2), (2.3) as the starting point. In order to simplify it further, we will 
switch to dimensionless variables (denoted with primes): 

, X . 77 , , h . U . t . C 

x=-, T] =-, h=-, u=—, t=— , andC=-, 
t a a -7=^ -£= a 

vgd \fgd 

where a and £ are the characteristic wave amplitude and wavelength correspondingly. We 
can compose three important dimensionless numbers which characterize the Boussinesq 
regime: 

e:= ~ « 1, ^ := (^) 2 « 1, and S := j 2 = 0(1), 

where S is the so-called Stokes-Ursell number [59], which measures the relative importance 
of dispersive and nonlinear effects. In unsealed variables the Peregrine- Wu system takes 
the following form (for simplicity we drop out the primes): 

r} t + ((h + eri)u) x = -h t , (2.4) 

up up- 

u t + euu x + r] x = —h(h t + (hu) x ) xt -—h 2 u xxt . (2.5) 

To simplify the problem, we will reduce the Boussinesq system (2.2), (2.3) to the unidi- 
rectional wave propagation. For instance, in the original work of T. Wu [62] a similar 
reduction to a forced KdV is also performed. However, the resulting model in our work 
will be of the BBM-type [3], since it possesses better numerical stability properties. 

The reduction to the BBM equation can be done in the following way [35]. The horizontal 
velocity u can be approximatively represented in unsealed variables as 

u = +i] + eP + fi 2 Q + 0(e 2 + efj 2 + fi A ), (2.6) 

where P(x,t) and Q(x,t) are some functions to be determined. The sign + in front of 
77 means that we consider the waves moving in the rightward direction. Substituting the 
representation (2.6) into unsealed Boussinesq equations (2.4), (2.5) yields two equivalent 
relations: 

7] t + j] x + eP x + \x 2 Q x + 2er]r] x - e((r]) x = 0(e 2 + e/2 2 + //), (2.7) 

up- up 

T] t + rj x + eP t + ii 2 Q t + er]r] x = —h{-C t + {hrj) x ) xt -—h 2 r] xxt + 0{e 2 + e^ 2 Mfi.' i i 
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By subtracting two last asymptotic relations we obtain the following compatibility condi- 
tion: 

e (P x - P t ) +f i 2 (Q x - Q t )+e mx -e (( V ) x = + n x ) xt+ ^h 2 r 1xxt+ 0{e 2 + e^+^). (2.9) 



For the right-going waves we have the following identities: 

P t = -P x + 0(e), Q t = -Q x + 0(e). 

Finally, the unknown functions P x and Q x can be chosen to satisfy asymptotically the 
compatibility condition (2.9), which yields 

2P X = ((r]) x - r]T] x , 2Q X = ~£ xtt - -rj xxt , 

where we used also the analytical representation of h(x, t) = l~e(^(x, t). The BBM equation 
in unsealed variables can be now easily obtained by substituting expressions for P x and Q x 
into equation (2.7): 

e u? (I 1 \ 

Vt + Vx + - (((v)x - VVx) + — [-Cxtt - -Vxxtj + 2erjr] x - s((r]) x = 0. 

Turning back to physical variables, the generalized forced BBM (gBBM) equation takes 
the following form: 

vt + ( y/7dn + yfl (lv 2 - \cv)) x - jvxxt - -\^ d Cxtt. 

In subsequent sections we will use this equation to model wave-bottom interaction. How- 
ever, in order to simplify the notations, we will introduce a new set of dimensionless 
variables, where all the lengthes are unsealed with the water depth d, velocities with \fgd 

fd 

and the time variable with \ l —. In this scaling the gBBM equation reads: 

V 9 

Vt + (v + -£) 2 ~ ^Cvj -^Vxxt = ~(xtt- (2.10) 

Recall that here r\ is the unknown free surface elevation, and ( is a given function, which 
is the topography of moving body defined by (2.1). The last equation has to be completed 
by appropriate initial and boundary conditions (when posed on a finite or semi-infinite 
domain) : 

rj(x,0) =?7o(x), xeR. (2.11) 

3. Well-posedness of the gBBM equation 

In this section we give a proof of the well-posedness of gBBM equation (2.10) in the 
Sobolev spaces H s := H S (R). First, we have the following result of the well-posedness of 
this problem. 

Theorem 1. For any ( e C 2 ([0, oo), H s ) nC([0, oo), H^ +1 ) and r] e H s , s > the problem 
(2.10), (2.11) admits a unique solution 77 6 C([0,oo),H s ). 
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Proof. Uniqueness. Let us assume that for some given functions C and 770 our problem 
(2.10), (2.11) admits two different solutions 771 and 77 2 . The difference 77 := 771 — 772 satisfies 
the following initial-value problem: 

Vt + (rj + -77(771 + 7/2) - -C?y) -g7fet = 0, 77(2, 0) = 0. (3.1) 

Taking the scalar product of (3.1) with 77 in L 2 , we obtain: 

M 

The Sobolev and Holder inequalities imply 



<C(||»7i||o + ||^||o+||C|o)||»7||i. 



Thus, the application of the Gronwall inequality yields 77 = 0. □ 
Existence. For any fixed time horizon T > let us show that our problem (2.10), (2.11) 
has a solution 77 e C([0, T), H s ). J. Bona & N. Tzvetkov (2009) [7] proved that for any 
given initial data 770 e H s the following BBM equation 

u t + u x + uu x - u xxt = 0, u(x,0) = u (x), 

admits a unique solution u e C([0, oo),H s ) (for s < they proved that the system is ill- 
posed). Using a scaling argument, we have also the well-posedness of the same equation 
with some positive coefficients: 

u t + {u + ^u 2 ^ --u xxt = 0, u(x,0) = rjo(x). (3.2) 

We seek a solution of (2.10), (2.11) in the form 77 = u + v, where u e C([0, 00), H s ) is the 
solution of (3.2) and v satisfies 

v t + (v + ^(v 2 + 2uv) - ^((v + - ^v xxt = ~Cxtt, v(x,0) = 0. 

Let us prove the existence of such v e C([0, 00), H s ) by induction on [s]. First, we 
assume [s] = 0. Taking the scalar product of (3) with v in L 2 , we obtain 

— f(-v 2 (x,t) + — vl(x,t)\ dx = - f —CxttV dx + f (-uv — C(v + u)\ v x dx. 
dtJ\2 K,J 12 xK ' ') J 4 s Jr\2 2 sv 7 

R K 

Using the Sobolev and Holder inequalities, we get 

* c(\\Cu(-,t)U\v(-,t)\\i + R,*)lloK-,*)ll? 

+ ||C(-,t)||iK-,i)||i(K-,i)||o + lk(-,t)||o)l (3.3) 
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After integrating (3.3) on the interval (0,t) we obtain 



T 

|2 



\\v(;t)\\i <C SUp \\v(;t)\\l / \\\Ctt(;s)\\o + H;s)\\oH; S )\\ 1 + 

IIC(-,*)lli (ll«(-,*)||o + ||«(-,*)||o) W (3.4) 

which is valid for any t e [0,T]. Hence, sup also can be estimated by the right 

u[o,T] 

hand-side of (3.4). Dividing by sup and applying the Gronwall inequality, we 

te[0,T] 

deduce 



sup \\v(;t)\\i<C [ (\\( tt (;s)\\a + K(;s)\\ 1 \\u(;s)\\ )ds> 
= rnn J 



te[0,T] 



exp^ (K-, S )|| + lC(-, S )|| 1 )d S 

Using this estimation and applying Picard iteration scheme, we obtain existence of the 
solution in the case s e [0, 1]. 

By induction, now we assume the existence of v for [s] < a for some integer a > 1 and let 
us prove it for [s] = a. To this end, let us take the 4-^, a < [s] derivative of (3), multiply 
the resulting equation by v a := 4-jjf and integrating in x over E, we obtain 

^ / (\vl(x, t) + L»£{x, t))dx = -f[\ (^U) v a ) dx + 

R R 

/ {l^ 2 + - \{% {C{V + U)) ) ^ (3 ' 5) 

R 

All the terms can be treated as above, except the term / ^;(v 2 )v ax dx, which is not zero 

R X 

in general. Using the induction hypothesis and the fact that a - 1 > 1, we can estimate 



Idx~> 2)v - 



dx 



< CWVaWllvW^ < M\\v a \\. 



Using the last estimation along with (3.5), the Sobolev and Holder inequalities, as above, 
we prove the required estimation for v a , which completes the proof. □ 

3.1. Optimization problem. In this section we turn to the optimization problem for the 
gBBM equation (2.10). We assume that the wavemaking piston is a solid, non-deformable 
object. Thus, its shape, given by a localized function Co(^), is preserved during the motion 
and it is sufficient to prescribe the trajectory of its barycenter only x = xo(t). Consequently, 
the time- dependent bathymetry is given by the following equation 

h(x,t) = d- Co{x - Xo(t)). 
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The piston shape Co( x ) an d its trajectory xo(t) will be determined as a solution of the 
optimization problem. More precisely, in the next section we will find numerically these 
functions in order to produce the largest possible wave (in L2 sense) in a given subinterval 
X = [a, b] of the numerical wave tank at some fixed time T > 0. In other words, we minimize 
the following functional: 

J(xoXo) = ~J~ r/(x,T) 2 dx — ► min, (3.6) 

where rj(x,t) is the solution of (2.10), (2.11). The existence of this solution is proven in 
the following 

Theorem 2. For any constants e, M > 0, there exists (xq,Q) 6 such that 

where Bm is a closed ball in H 2+£ [0,T] x Hq +e ([0, 1]) centered at origin with radius M . 

Proof. Let (£q , Co) t> e an arbitrary minimizing sequence of J. Since H 2+£ [0, T]xHq +£ (\0, 1]) 
is Hilbert space, extracting a subsequence, if it is necessary, we can assume that there is 
(4, Co) e B M such that (a#,C?) (xq, Q) in B M . 

Let us denote rf 1 the solution of (2.10), (2.11) with C = C n := CoC^ ~ x o (0)- Let us show 
that we have r/ n (T) -*■ rj*(T) in L 2 , where r]* is the solution of (2.10), (2.11) with C = C*- 
Indeed, for fj n ' m := rf 1 - r] m we have 

(3.7) 

Since C t " -> C« := d u(Q( x ~ x o( f ))) in -^ 2 ([0, T], L 2 ), multiplying (3.7) in L 2 by ^> m , 
integrating by parts and applying the Gronwall inequality, we obtain that rf 1 is a Cauchy 
sequence in if 1 . Hence, 

J(x ,C *) = lim J(*o\C n )- 

n->oo 

This completes the proof of the theorem. □ 

However, in practice, the functional (3.6) has to be completed by appropriate constraints 
in order to provide a solution realizable in practice. For example, the speed of the under- 
water piston is limited by technological and energy consumption constraints. Some more 
realistic formulations will be addressed numerically in the next Section. 

4. Numerical results 

In order to discretize the gBBM equation (2.10), posed on a finite interval we use 

a modern high-order finite volume scheme with the FVCF flux [29] and the UN02 recon- 
struction [33]. The combination of these numerical ingredients has been extensively tested 
and validated in the context of the unidirectional wave models [22] and Boussinesq-type 
equations [20, 21]. For the time-discretization, we use the third-order Runge-Kutta scheme, 



~n,m , / ~n m . 

Vt + 1 V + 
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Parameter Value 

Computational domain [a, /3]: [-5, 10] m 

Wave quality evaluation area [a, b\. [0, 6] m 

Number of discretization points N: 1000 

CFL number: 1.95 

Gravity acceleration g: 9.8 ms~ 2 

Undisturbed water depth d: 1.0 m 

Final simulation time T: 7.0 s 

Piston motion total time Tf. 3.0 s 

Piston length £ : 1.0 m 

Piston maximal height ao/d: 0.1 

Piston starting point x = 0.0 m 
Upper bound of the piston position x max : x max = 4.5 m 

Upper bound of the piston speed vf. 1.5 m/s 

Wave generation limit Xf\ x = 1.0 m 

AA-wave solution ceter x m : x = 2.0 m 



Table 1. Values of various parameters used in numerical computations. 



which is also used in the ode23 function in Matlab [52]. In all experiments presented below 
we assume that the water layer is initially at rest: 

t](x,0) = r} (x) = 0. 

The computational domain [a,/3] is discretized in N equal subintervals, called usually the 
control volumes. The time step is chosen locally in order to satisfy the following CFL 
condition [14] used in shallow-water models: 



max Ui + \fgd 

\<i<N 

The values of all physical and numerical parameters used in simulations are given in Table 
1. 

On the left and right boundaries we apply the Neumann-type boundary conditions which 
do not produce reflections. In any case, we stop the simulation before the generated wave 
reaches the right boundary. We recall that the gBBM equations (2.10) describes the 
unidirectional (rightwards, for instance) wave propagation. So, the influence of the left 
boundary condition is negligible. 

Let us describe the constraints that we impose on the shape (,§{x) > and trajectory 
xo(t) of the underwater wavemaker. First of all, we fix the length 2£ Q of this object. Then, 
we assume that its height is also bounded: 

Co(x) ^ 
max — - — < ao • 

xeR d 
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We allow the piston to move during the first Tf s, its motion always starts at the same 
initial point x% and it is confined to some wave generation area [xo(0),Xf] £ [a, f3\. 

support) £ [0,7)], x (0) = x° , x° <x (t)<x f , V*e[0,7>]. 

However, the cost function J(x ,Co) is evaluated at time T > Tf so that the generated 
waveform can evolve further into the desired shape. 

Moreover, we require that the piston speed and acceleration are bounded, since too 
fast motions are difficult to realize in practice because of the gradually increasing energy 
consumption: 

sup (\x' (t)\ + J-\x%(t)\)<v f 

U[0,T f ] V 9 

In order to parametrize the wavemaker shape, we use only three degrees of freedom 
£ ( - lf)' Co(0), Cod 1 ) which represent the height of the object in three points equally 
spaced on the suppCo- Finally, the continuous shape is reconstructed by applying the 
interpolation with cubic splines 1 through the following points: 

Ho,o), (-|,Co(-|)), (o,Co(o)), (|,Co(|)), (4,0). 

In a similar way we proceed with the parametrization of the piston trajectory x (t) which 
is represented with 4 degrees of freedom (three in the interior of the interval (0,7)) and 
the final point x (Tf) which is not fixed as in the case of Co)- Obviously, more degrees of 
freedom can be taken into account when it is needed for a speecific application. However, 
the number of degrees of freedom determines the dimension of the phase space where we 
find the optimal solution. In examples below we operate in a closed subset of M7. In order 
to obtain an approximate solution to our constrained optimization problem, we use the 
function f mincon function of the Matlab Optimization Toolbox. This solver is a gradient- 
based optimization procedure which uses the SQP algorithm. The iterative process is 
stopped when the default tolerances are met or the maximal number of iterations is reached. 
The values of default tolerances are TolCon (tolerance on the constraint violation) = TolFun 
(tolerance on the function value)=TolX (tolerance on the parameter estimates)=l(T 6 and 
Maxlter (maximum number of iterations) =200. 

As the first numerical example, we minimize the functional J(xq,Co) subject to con- 
straints described in above. Basically, this cost-function measures the wave deviation from 
the still water level in a fixed portion [a, b] of the wave tank, consequently, bigger waves 
in this interval will provide lower values to the functional J. The result of the numerical 
optimisation procedure is represented on Figure 3. The free surface elevation computed 
at the final time T is shown on Figure 3(a). One can see that in the region of interest 
[a, b] = [2,4] m we have a big depression wave which is followed by a wave of elevation. We 
make a conclusion that we succeeded to generate a wave suitable for surfing purposes in 
artificial environments. The computed shape of the underwater object is shown on Figure 
3(b) and its trajectory is represented on Figure 3(c). It is interesting to note that the com- 
puted optimal shape is composed of two bumps. The piston trajectory can be conditionally 

1 Cubic splines ensure that the interpolant belongs to the class C 2 . 
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Figure 3. Computed numerically the optimal piston shape and its trajec- 
tory which minimize the functional J(xq,Co)- 



decomposed into three parts. During the first 1.25 s we have a stage of slow motion, which 
is followed by a rapid acceleration and, during the last 0.75 s, we can observe a slight 
backward motion of the piston before it is frozen in its final point. Then, the wait has 
T - Tj = 4 s to evolve before its quality is estimated using the cost function J(x , Co)- 

Since the choice of the functional to minimize is far from being unique, we decided to 
perform some additional tests. Instead of maximizing the wave height, one can try to 
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maximize, for example, the wave steepness in a given portion of the wave tank. In other 
words, we will minimize the following functional (subject to the same constraints as above): 

Ji(xo,Co) = -f Vx(x,T)dx. 

i 

The result of the numerical optimization procedure is shown on Figure 4. One can see on 
the free surface snapshot 4(a) that effectively the wave became steeper. The optimal shape 
of the wavemaker is almost the same as for the functional J(xo,(o)- However, the piston 
trajectory is almost monotonic and close to the uniform motion. This solution might be 
easier to implement in practice. 

Finally, we can simply minimize the mismatch between the obtained solution and a fixed 
desired wave profile: 

M x o, Co) = J (r](x,T) -r] T (x)) 2 dx, 
i 

where t]t(x) is a given function on the interval X. To illustrate this concept, in numerical 
computations we take the AA-wave ansatz put forward by S. Tadepalli & C. Synolakis 
(1994, 1996) [53, 54]: 

The first profile T]^ (x) corresponds to the leading elevation AA-wave solution (LEN), while 
the second function rjj, (x) is a typical leading depression N -wave (LDN). The results of 
optimization procedures are shown on Figures 5 and 6. One can notice that the resulting 
optimal shapes of the wavemaker are completely different (see Figures 5(b) and 6(b)). For 
the surfing applications the LDN wave might be more interesting. It requires also more 
uniform piston motion comparing to the LEN wave (see Figures 5(c) and 6(c)). 

Remark 1. The arguments used to prove Theorem 2 can be also applied to show the 
existence of minimizers for the functionals Ji(xo,Co) an( ^ ^( x oXo)- 



5. Conclusions 

In the present work we considered the water wave generation problem by disturbances 
moving along the bottom. This problem has many important applications going even to 
the design of artificial surfing facilities [34]. In order to study the formation of water waves 
due to the motion of the underwater piston, we derived a generalized forced BBM (gBBM) 
equation. The existence and uniqueness of its solutions were rigorously established. The 
trajectory of the piston is determined as the solution of a thoroughly formulated opti- 
mization problem. The existence of minimizers is also proven. Finally, the theoretical 
developments of this study are illustrated with numerical examples where we solve several 
constrained optimization problems with various forms of the cost functional. The resulting 
solutions are compared and discussed. 

In future studies this problem will be addressed in the context of more complete bidi- 
rectional wave propagation models of Boussinesq-type [5, 17, 39, 16]. The optimization 
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Figure 4. Computed numerically the optimal piston shape and its trajec- 
tory which minimize the functional Ji(xq,Co)- 



algorithm can be also further improved by evaluating the gradients analytically, for exam- 
ple. From physical point of view, one may want to include some weak dissipative effects 
for more realistic wave description [18]. 
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Figure 5. Computed numerically the optimal piston shape and its tra- 
jectory which minimize the functional ^2(^0, Co) and the terminal state 
Vt\ x ) = (x-x m )sech 2 (x-x m ). 
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